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We present a measurement of the CP-violating asymmetry in decays using 535 million 

BB pairs collected with the Belle detector at the KEKB e+e^ collider. We measure CP-violating 
coefficients A = 0.16 ± 0.21 (stat) ± 0.08 (syst) and S = 0.19 ± 0.30 (stat) ±0.08(syst). These 
values are used to determine the unitarity triangle angle (j>2 using an isospin analysis; the solution 
consistent with the Standard Model lies in the range 54° <(f>2< 113° at the 90% confidence level. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



CP violation in the Standard Model is attributed 
to the presence of an irreducible complex phase in the 
Cabibbo-Kobayashi-Maskawa (H (CKM) quark -mixms; 
matrix. The unitarity of the CKM matrix leads to six 
triangles in the complex plane. One such triangle is 
given by the following relation among the matrix ele- 
ments: VudV:t, + VcdV:t, + VtdV*t, = 0. The phase angle </.2, 
defined as arg[— (VtdV^l)/ (KjrfTCb)] ' '^^^ t)e determined by 
measuring a time-dependent CP asymmetry in 6 — > uud 
decays such as B'^ tt+tt^, (/Ott)", and p^p~ The 
time-dependent rate for B — > P~ decays tagged with 
B^{q = -1-1) and B^ {q = —1) mesons is given by 

g-|At|/rgo 

Vpp{At) — {1 + g[^cos(AmAt) (1) 

4tbo 

+Ssm{AmAt)]}, 

where Tgo is the B'^ lifetime, Am is the mass difference 
between the two B° mass eigenstates, At is the proper- 
time difference between the two B decays in the event, 
and A and S are CP asymmetry coefficients. If the decay 
amplitude is a pure CP-even state and is dominated by a 
tree diagram, S = sin(2(/)2) and ^ = 0. The presence of 
an amplitude with a different weak phase (such as from a 
loop diagram) gives rise to direct CP violation and shifts 
S from sin(202)- However, the size of a loop amplitude is 
constrained to be small by the small branching fraction 

of PO^PV i- 

The CP-violating parameters receive contributions 
from a longitudinally polarized state (CP-even) and two 
transversely polarized states (an admixture of CP-even 
and CP-odd states). Recent measurements of the po- 
larization fraction by Belle Q and BaBar Q show that 



the longitudinal polarization fraction is near unity (/^ = 
0.968 ±0.023 [i]). 

Here we present an improved measurement of the CP- 
violating coefficients A and S using 492 fb~^ of data con- 
taining 535 million BB pairs. This data sample is about 
a factor of two larger than that used in our earlier publi- 
cation 0]. In addition we have modified the event selec- 
tion by reducing the threshold on a continuum suppres- 
sion variable; this increases our reconstruction efficiency 
by about 70%. We subsequently introduce a probability 
density function (PDF) for the continuum suppression 
variable into the likelihood function; this provides ad- 
ditional discrimination between signal and backgrounds. 
The expected improvement in the statistical error of A 
and S due to the new event selection is about 12%. 

The BB pairs were collected with the Belle detector 0] 
at the KEKB Q e+e" asymmetric-energy (3.5 GeV on 
8.0 GeV) collider with a center-of-mass (CM) energy at 
the T(4S') resonance. The T(4S') is produced with a 
Lorcntz boost of = 0.425 nearly along the z axis, 
which is oriented antiparallel to the positron beam. Since 
the P° and P° mesons are produced approximately at 
rest in the T{4S) CM system, the decay time difference 
At is related to the distance between the decay vertices 
of the two B mesons as At ~ Az/Pjc, where c is the 
speed of hght. 

The Belle detector [3] is a large-solid-angle spectrome- 
ter. It includes a siHcon vertex detector (SVD), a 50-layer 
central drift chamber (CDC), an array of aerogel thresh- 
old Cherenkov counters (ACC), time-of-flight scintilla- 
tion counters (TOF), and an electromagnetic calorime- 
ter (ECL) comprised of CsI(Tl) crystals located inside a 
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FIG. 1: Left: projections in M^c for events satisfying 
-0.10 GeV < A_E < 0.06 GeV. Right: projections in AE 
for events with 5.27 GeV/c^ < A/bc < 5.29 GeV/c^ The top 
plots correspond to good quality tags (0.75 < r < 1.0), and 
the bottom plots correspond to lower quality tags (r<0.75). 
The curves show fit projections: dashed is p'^ + ptvtt, dot- 
ted is qq, dot-dashed is b—^c, small solid is b-^u, and large 
solid is the total. For these plots the TZ requirement has been 
tightened to increase the ratio of signal to background. 



superconducting solenoid coil that provides a 1.5 T mag- 
netic field. 

We reconstruct p+p^ decays by combining two 

oppositely charged pion tracks with two neutral pions. 
Each charged track is required to have a transverse mo- 
mentum prp > 0.10 GeV/c in the laboratory frame and 
originate within dr < 0.2 cm in the radial direction and 
within |o?z| <4.0 cm along the z-axis from the interaction 
point, which is determined run- by-run. A track is iden- 
tified as a pion using information from the CDC, ACC 
and TOF systems. Tracks matched with clusters in the 
ECL that are consistent with an electron hypothesis are 
rejected. 

The 7r° candidates are reconstructed from 77 pairs 
with an invariant mass in the range 117.8 MeV/c^ < 
M^^ < 150.2 MeV/c^ (about ±3a in m^o resolution). 
Photons are required to have energy > 50 MeV in the 
ECL barrel region (32° <e< 129°) and > 90 MeV in 
the endcap regions (17° < 6* < 32° and 129° < 6* < 150°), 
where 9 denotes the polar angle with respect to the z 
axis. 

To reconstruct mesons, we combine tt^ candidates 



with TT candidates. 



The TT^TT*^ combination must have 



an invariant mass in the range 0.62 GeV/c^ < A/^i^rO < 
0.92 GcV/c^. To reduce combinatorial background, we 
reject tt^'s with p < 0.35 GcV/c in the CM frame. We 




FIG. 2: TZ distribution for high-purity tagged events satisfying 
5.27 GeV/c^ < Mbc < 5.29 GeV/c^ and -0.10 GeV < AE < 
0.06 GeV. The curves show fit projections: dashed is p^p~-f 
pnn, dotted is qq, dot-dashed is b ^ c, small solid is b~* u, 
and large solid is the total. 



also require — 0.80 < cos 6|£ < 0.98, where 9^. is the angle 
between the direction of the 7r° from the and the 
negative of the momentum in the p^ rest frame. 

B*^ p^ p^ decays are identified using the beam- 
energy-constrained mass M^^ 



beam 



Pb 



ergy difference AE' = E^ 



^beam' ^hcrC E^^^^^ 



is the 



beam energy, and E^ and pb are the energy and mo- 
mentum of the reconstructed B candidate, all evaluated 
in the CM frame. 

The flavor of the B meson accompanying the i?" 
p'^ p~ candidate is identified via its decay products: 
charged leptons, kaons, and A's. A tagging algorithm Q 
yields the flavor of the tagged meson, g, and a quality 
factor, r. The parameter r ranges from for no flavor 
discrimination to 1 for unambiguous flavor assignment. 
We divide the data sample into six r intervals (denoted 
^=1,2,. ..,6). The wrong tag fractions cu^ for these inter- 
vals and the differences Aw^ in these fractions between 
B'^ and B^ decays are determined from data 

The dominant background originates from e~^e~ — > 
qq {q = u,d,s,c) continuum events. To separate qq jet- 
like events from more spherical BB events, we use event- 
shape variables, specifically. 16 modified Fox-Wolfram 
moments 11 1 combined into a Fisher discriminant [l^ . 
We form signal and background likelihood functions 
and CgQ by multiplying the PDF for the Fisher discrim- 
inant by a PDF for cos6'^, where 9^ is the polar angle in 
the CM frame between the B direction and the beam 
axis. The PDFs for signal and qq are obtained from 
Monte Carlo (MC) simulations and the data sideband 
5.23 GeV/c2<Afbc<5.26 GeV/c^, respectively. We cal- 
culate the ratio TZ = Cg/{Cg + C^q) and make a loose 
requirement 7^ > 0.15. 

The decay vertices of a p^ p^ candidate and the tag- 
side B meson are reconstructed using charged tracks that 
have a sufficient number of SVD hits and an interaction 
point constraint. The vertex reconstruction algorithm is 
described in Ref. [l^]. 

The analysis is organized in two main steps. We first 
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determine the yields of signal and background compo- 
nents from a fit to the three-dimensional {M^^, AE,TZ) 
distribution. Here, candidates are required to satisfy 
5.23 GeV/c2 < M^^^ < 5.29 GeV/c^, -0.2 GeV < AE < 
0.26 GeV, and TZ > 0.15. We subsequently perform a 
fit to the At distribution to determine the CP param- 
eters A and <S. The signal region used for the At fit is 
5.27 GeV/c2 < M^^ < 5.29 GeV/c^, -0.12 GeV < AE < 
0.08 GeV, and 7e>0.15. 

About 12.6% of events contain multiple B° — > p'^p^ 
candidates, most of which arise from fake tt'^'s combining 
with good tracks. We select the best candidate based on 



the TT masses, i.e., minimizing "^^o (r 



TOtto) . For 



the small fraction (3%) of multiple-candidate events that 
arise due to extraneous tt* tracks combining with a single 
TT*^, wc select one randomly. Signal decays that have at 
least one tt meson incorrectly identified are referred to as 
"self-cross-feed" (SCF) events. 

The likelihood function used to determine the event 
yields is given by 



iVe, 



N,Vj{Ml,,AE\n^) 



,(2) 



where j indicates one of the following event categories: 
signal and /otttt non-resonant decays, SCF events, contin- 
uum background [qq), b ^ c background, and charmless 
{b — s- u) background. Nj is the yield of each category, 
Vj(Ml^, AE\ TV) is the PDF for the «-th event for cate- 
gory j , and A'ovt is the total number of events in the fit. 
Except for the small contributions of 6 — > w background 
and SCF events, the yields Nj are determined from the 
fit. Due to the similar shapes of the Mj^^, AE, and TZ 
distributions for signal and pirir events, we cannot distin- 
guish these two components; the fraction of p-n-n events is 
measured in Ref. Q and constitutes (6.3 ± 6.7) % of the 
total Np+p~^p^^ signal. The fraction of SCF events is de- 
termined from MC simulation. The M^^^ and AE shapes 
for the signal and SCF components are modeled by a two- 
dimensional smoothed histogram obtained from a large 
MC sample. To take into account a small difference be- 
tween the MC and data, the M^^ - AE shapes are cor- 
rected according to calibration factors determined from 
a S+ ^ D°p+, Z)o ^ A'+tt-tt" control sample. The 7^ 
shapes are modeled by one-dimensional histograms, also 
obtained from MC simulation. 

The PDF for 6 — > c background is the product of a 
threshold ARGUS function [31 for M^^, a quadratic 
polynomial for AE, and the sum of a Gaussian and a 
third-order polynomial for TZ. The shapes of the AE 
and TZ distributions depend on the tag quality bin £. Pa- 
rameters for all distributions are obtained from the MC. 

The and AE PDFs for qq are modeled by an AR- 
GUS function and a linear function, respectively. The 
AE slope depends on TZ and the tag quality bin £. The 



shape parameters for Mj^^ and AE are floated in the fit. 
The TZ PDF for qq background is taken to be an eighth- 
order polynomial function; the coefficients depend on the 
bin £ and are determined from a data sample collected at 
a CM energy - 60 MeV below the T(4S'). 

The b ^ u background is dominated by _B ^ 
(pTT, a^TT, a^p) decays. We estimate the i?^ — > (ain)^ 
branching fractions (which are unmeasured) to be (20 ± 
10) X 10~^ using the measured value for B^^af tt^ [l^ . 
For S± ^ (aip)^ we assume branching fractions of 
(30 ± 15) X 10^^, consistent with the present upper limit 
for SO^a^pT (<6 X 10"^ [H). The fraction of b-^u 
events is very small (0.37%) and thus is fixed in the fit 
according to the prediction of MC simulation. A fit to 



176843 events maximizing £ yields N, 



PP+pTTTT 



= 576 ± 53. 



Figures [T] and [5] show the M^^, AE, and TZ distributions 
along with projections of the fit result. 

The CP-violating parameters A and S are obtained 
using an unbinned ML fit to the At distribution. The 
likelihood function for event i is given by 

A=E / fni^^)KWR,,iAt\At')dAt', (3) 



where n is one of the six event categories: correctly 
reconstructed signal, SCF events, pmr non-resonant 
events, b c background, qq background, and b 
u background. The weights /„ are functions oi x £ 
{M^^, AE,TZ) and are normalized to the event fractions 
obtained from the {M^^,AE,TZ) fit. The PDFs 'P„(At) 
are convolved with the corresponding At resolution func- 
tions i?„. Both /„ and ^^(Ai) depend on the tag quality 
bin £. 

The signal PDF is given by Eq. (1) modi- 
fied to take into account the effect of incor- 



rect flavor assignment: 



|At|/r^o 



/(4. 



{1 - qAujf + q{l - 2ijjf) [Aco8{Am At) + 5 sin(Am Ai) ]}. 
As the fraction of longitudinal polarization is close 
to 100%, we assume that A~Al, S = Sl, and consider 
the potential contribution from a transversely polarized 
amplitude as a systematic uncertainty. The signal PDF 
is convolved with the same At resolution function as 
that used for Belle's sin 20]^ measurement [loj . 

The fraction of SCF events with incorrectly recon- 
structed vertices is estimated from MC simulation to be 
(6.5 ± 0.1)% of aU signal events. The PDFs Vp^^ and 
T^SCF exponential with T = Tg and t«0.96 ps (from 
MC), respectively; these are smeared by a common reso- 
lution function. 

The At PDFs for the backgrounds are modeled as 
a sum of prompt and exponential components: T'k = 
fgS{At) + (1 - /|)e^l^*l/'^''/2rfc, where k represents con- 
tinuum, b ^ c, and b ^ u backgrounds, /| is the fraction 
of the prompt component, S{At) is the Dirac delta func- 
tion, and Tfe is an effective lifetime. These PDFs are con- 
volved with a resolution-like function Rk parameterized 
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FIG. 3: The At distribution and projections of the fit for 
events satisfying 0.5 < r < 1.0: (a) g = +1 tags, (b) g = — 1 
tags. The hatched region shows signal events. The raw CP 
asymmetry is shown in (c). For these plots the 71 require- 
ment has been tightened to increase the ratio of signal to 
background. 
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FIG. 4: 1 — C.L. vs. (j)2- The horizontal lines denote C.L. = 
68.3% (solid) and C.L. = 90% (dashed). 



error due to transverse polarization is obtained by first 
setting f]^ equal to its central value and varying Aj^, 
Srp from —1 to +1; then, conservatively assuming that 
the transversely polarized component (with fraction /t = 
1 — Jl) is pure CP-odd for which At = Al^ St = —Sl, 
and varying by its error. Summing up in quadrature 
all systematic uncertainties, we obtain overall systematic 
errors of ± 0.08 for both A and S. Thus, 



as a sum of two Gaussian functions. Parameters for Vk 
and Rk are determined from a data sideband for contin- 
uum background and from large MC samples for b ^ c 
and b u backgrounds. To account for small correla- 
tions between the shape of the At distribution and TZ for 
qq background, the parameters are obtained separately 
for low (0.15 < 7^ < 0.75) and high (0.75 < 7^ < 1.0) 7^ 
regions. 

We determine A and S by maximizing log , where 
i runs over the 18016 events in the {M^^, AE,TZ) sig- 
nal region. The results are ^ = 0.16 ± 0.21 and S = 
0.19 ± 0.30, where the errors are statistical. The corre- 
lation coefficient is —0.10. These values are consistent 
with no CP violation {A — S — 0); the errors are consis- 
tent with MC expectations. Figure [3] shows the data and 
projections of the fit result. 

The sources of systematic error are listed in Table HI 
The error for most sources is evaluated by varying the 
corresponding parameters by ±1 standard deviation (a). 
The effect of a possible asymmetry in 6 — > c and qq is 
evaluated by adding such an asymmetry to the 6— >c and 
qq At distributions. The uncertainty due to a possible 
asymmetry in pmr non-resonant decays is estimated by 
varying Ap-KTr and vSp^Tr by 0.68, corresponding to a 68% 
confidence interval of a free distribution. We vary the 
branching fractions for aip and aitt decays and also allow 
for a CP asymmetry of up to 100% in these modes. The 



= 0.16 ±0.21 (stat) ±0.08 (syst) (4) 
S = 0.19 ±0.30 (stat) ±0.08 (syst). (5) 

These values are consistent with, and supersede, our pre- 
vious measurement They are also consistent with 
results obtained by BaBar 

We constrain using an isospin analysis [l7j . which 
allows one to relate six observables to six underlying pa- 
rameters: five decay amplitudes for B ^ pp and the 
angle 02. The observables are the branching fractions 
for B p+p-, p+p^ 0, and S; the CP param- 
eters A and iS (our results); and the parameter ApOpO 
for B p^ p^ decays. The last parameter is not yet 
measured, but nevertheless one can constrain 02- The 
branching fractions must be multiplied by the corre- 
sponding longitudinal polarization fractions 0]. We ne- 
glect possible contributions from electroweak penguins 
and 7=1 amplitudes [31 and possible interference be- 
tween signal and non-resonant components. We follow 
the statistical method of Ref. [l^ and construct a ^{^-i) 
using the measured values and obtain a minimum (de- 
noted Xmin); we then scan (/jj from 0° to 180°, calculating 
the difference Ax^ = x^(02) ^ Xmin- insert A^ into 
the cumulative distribution function for the x^ distribu- 
tion for one degree of freedom to obtain a confidence level 
(C.L.) for each value. The resulting function 1 — C.L. 
(Fig. HI) has more than one peak due to ambiguities that 
arise when solving for (^^ ■ The "flat-top" regions in Fig. 2] 
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TABLE I: Systematic errors for CP coefficients A and S. 



Type 



(xlO-^) AS (xlO- 





+cr 


— a 


+cr 


— a 


Wrong tag fractions 


0.5 


0.5 


0.8 


0.8 


Parameters '^b^ 


0.2 


0.3 


0.6 


0.7 


Resolution function 


1.4 


1.5 


1.0 


1.7 


Background At distributions 


0.5 


0.5 


1.0 


1.1 


Component fractions 


1.5 


1.9 


3.9 


3.7 


pTvn non-resonant fractions 


1.2 


1.0 


1.5 


1.2 


SCF fraction, At PDF 


0.2 


0.2 


0.1 


0.1 


oiiape oi /V LUr 


n 8 

u.o 


n 7 


1 9 
r 


1 


Vertexing 


2.1 


2.1 


1.0 


1.3 


Possible fitting bias 


0.3 


0.0 


0.3 


0.0 


Background asymmetry 


1.1 


0.0 


0.0 


0.4 


b u asymmetry 


2.4 


2.9 


2.4 


3.2 


pTTTT asymmetry 


4.6 


4.6 


4.6 


4.6 


Transverse polarization 


3.8 


2.8 


4.6 


2.7 


Tag-side interference [16] 


3.7 


3.7 


0.1 


0.1 


Total 


8.3 


8.0 


8.4 


7.9 



arise because ApOpO is not measured. The solution con- 
sistent with the Standard Model is 62° < <j>2 < 106° at 
68% C.L. or 54° < < 113° at 90% C.L. Recently, an 
alternative model-dependent approach to extract (/)2 us- 
ing flavor SU{3) symmetry has been proposed [2^. This 
method could potentially give more stringent constraints 
on 02- 

In summary, we present an improved measurement of 
the CP-violating coefficients A and S in ^ p'^ p~ de- 
cays using 492 fb"-'^ of data, which corresponds to 535 
million BB pairs. These measurements are used to con- 
strain the angle (/)2. 
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